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Abstract 

Self-energy corrections of order O(o 2 ) to the parapositronium and or- 
thopositronium decay rates are calculated. Numerical values of the corre- 
sponding coefficients are B p = —3.74, B Q = 2.02. 
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Investigation of the parapositronium and orthopositronium decay rates is one 
of the important tasks in the QED bound state problem. Especially interesting 
experimental situation came into existence about orthopositronium decay width 
[|l|. ^|. |3|, U . While the experiments performed by the University of Michigan group 
[§, [J lead to strong disagreement between theoretical and experimental values of 
orthopositronium decay rate (6-9 standard deviations), in the experiments of Tokyo 
University group this discrepancy keep within the experimental accuracy f7|. The 
theoretical expression for the orthopositronium decay rate, known at present, can 
be written in the form: 
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A = -10.286606(10). 



Numerical value of (1) without a 2 correction is equal to 7.0382 /is 1 . Experimental 
results for the orthopositronium lifetimes are: 



T exp (o-Ps, gas measurement) = 7.0514(14) fis 1 [5], 

T exp (o — Ps, vacuum measurement) = 7.0482(16) fis^ 1 
T exp (o - Ps, Si0 2 measurement) = 7.0398(29) fis' 1 [7], 



(2) 



(3) 



The coefficient B Q should be equal to 41 |§ so that the expressions (1) and (3) 
coincide. Such value of B D can be derived on the basis of standard QED calculations 
[0, p. During some years the different QED effects of order 0(a 2 ) in the V th were 
studied [|H], [Tl|, [12|, [13|, [14|, [15] , but exact value of coefficient B Q is not known up to 
now. In the case of parapositronium the corrections a 2 were calculated recently in 
H], [17| . The theoretical expression for the decay rate of singlet positronium is |17| : 
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where B p = 1.75(30). The corresponding experimental quantity |T8 



T exp {p - Ps) = 7990.9 ± 1.7 us' 
agrees well with the value from (4): 

F th (p - Ps) = 7989.50^/s- 1 



(4) 
(5) 

(6) 



Two-loop contributions to the o-Ps lifetime may be devided into several classes: 
vertex, vacuum polarization, self-energy and annihilation corrections. Two-loop cor- 
rections due to vacuum polarization were calculated in [[□], |T2|. Some annihilation 
type contributions to V th (o — Ps) were studied in [14] , |i~5|1 . The second order electron 
self-energy corrections ([[!!], ^(J, ^2[ |23|) represent the special set of a 2 corrections 
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Pi p 2 p 3 

Figure 1: Feynman diagrams representing two-loop self-energy corrections to para- 
positronium decay rate. 



in the decay width of positronium. Such corrections can be studied independently by 
means of renormalized expression for the mass operator. In this work we calculate 
o? corrections to the positronium decay rate, connected with electron self-energy 
insertions to inner electron lines. We don't consider here self-energy vacuum po- 



larization corrections, which were calculated in [11], [12]. Corresponding Feynman 
diagrams are shown on Fig. 1-2. 

The decay width of o-Ps into three photons can be written in the form |9|, 2~4}| : 
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and amplitude M(o — Ps — > 37) has the following integral representation: 

d 3 p 



M(o -Ps^ 37) 
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where P=(M=2m,0) is the four-momentum of the orthopositronium, k\, k 2 , k 3 are 
photon momenta, ifi(p) is the positronium Coulomb wave function, M(ei,p) is the 
amplitude of the annihilation e + + e~ — > 37, II is the projection operator on the 3 5'i 
electron-positron state: 

is the polarization vector of orthopositronium. In the case of parapositronium we 
must change e — > 75. 

As it is well known the mass operator has no infrared divergences in the Fried- 
Yennie gauge (FY) |25| . So, we can perform the usual mass-shell subtraction without 
introduction of the photon mass. The bare mass operator has the following structure: 



E(p) = 5m + (l - Z 2 Y ) (p- m ) + S (R) (p). 



'IF 
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Then the self-energy correction to the electron line is determined by the substitution: 

1 1 , 1 



p — m p — m 



p — m 



(12) 



The renormalized self-energy operator of the second order in the FY gauge can be 
presented as follows [P 
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We have used this expression to obtain loop-after-loop contributions, presented by 
Pi and Oi graphs. Feynman diagrams P2, O3 contain one internal and one exter- 
nal loops. To find contributions of these diagrams to the mass operator we take 
expression (13) for inner loop and again use FY gauge for external photon. In- 
troducing Feynman parameterization, momentum cutoff and subtracting ultraviolet 
divergences in accordance with (11), we have obtain the following renormalized ex- 
pression for mass operator with two non-crossed photons: 
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Explicit expression for the overlapping two-loop self-energy diagram in the FY gauge 
was derived in [27|. Corresponding renormalized part of the mass operator has the 
form of integral representation over five Feynman parameters (functions a p , a m ): 
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We used this relation to calculate contributions of diagrams P3, O4. The results of 
our calculation of a 2 self-energy corrections to the p-Ps and o-Ps decay rates are 
shown in Tables 1,2. 

Table 1. Self-energy a 2 corrections to p-Ps decay rate. 
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Total contributions of the diagrams Pi 
are equal: 



P3; Oi — 4 to the coefficients B p , B Q 
-3.74, B S Q E = 2.02. (16) 
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Figure 2: Feynman diagrams representing two-loop self-energy corrections to or- 
thopositronium decay rate. 
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Bq E to the o-Ps lifetime increases the quantity 
' '"-~ [ Our results for finite part of self-energy 



The contribution of the term ^ 
7.0382 ps~ l by the value 8 • 1(T" fia 

contribution to the p-Ps decay rate do not coincide with the results of 



We 



think that the reason of such difference lies in the choice of gauge for exchanged 
photons and in the used renormalization scheme. The calculation of the modulus 
squared of the amplitude \M\ 2 was done by means of the system of analytical calcu- 
lations "Form" pSfl . Many self-energy corrections for Ps decay width, presented in 
the Tables were obtained by the numerical multidimensional integration by means 
of system " M athematica" | 29fl . Integral representations for the corrections to the 



orthopositronium coefficient B are given in the Supplement. 



Table 2. Self-energy a 2 corrections to o-Ps decay rate. 



Feynman diagram 
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A. Contributions of Feynman diagrams Oi — O4 to the co- 
efficient B. 

B °i = i7~ 2 — ST / dpi / dp2 Ta \ x 17 

8{n 2 -9)Jo J2-P1 Pip 2 (4 - pi - p 2 ) 

x [8 - 20 Pl + ~ P i(U Pl + 15p 2 ) - -pl(2 Pl + 5p 2 ) - -^Mpi + 2p 2 )], 

1 - Pi V 1 - Pi J 

27 f 2 A f 2 a (1-Pi + Pilnpi) 2 , 1B . 
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+\p\p2 + y/A + -^pI + ^p\p2 + \p\p\ + ^P 2 lP 2)+ 

T 33 X XX X5 X 
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